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Internal  solitary’  waves  induced  by  flow  over  a  ridge: 

With  applications  to  the  northern  South  China  Sea 

Ping-Tung  Shaw,1  Dong  Shan  Ko,2  and  Shenn-Yu  Chao1 
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[i]  The  generation  of  internal  solitary  waves  by  barotropic  tides  over  a  ndge  is  studied  in 
a  nonhydrostatic  numerical  model  under  idealized  oceanographic  settings.  The 
experiments  examine  the  effects  of  ndge  width,  barotropic  tidal  strength,  and  stratification 
on  wave  generation.  The  barotropic  tidal  flow  produces  internal  wave  beams  emitting 
from  the  ridge  top  if  the  slope  of  the  ridge  exceeds  a  critical  value  equal  to  the  slope  of  the 
wave  beam.  Reflection  and  refraction  of  a  wave  beam  in  an  upper  ocean  waveguide 
associated  with  a  strong  shallow  thermoeline  produce  honzontally  propagating  internal 
tides.  When  the  local  Froude  number  over  a  ridge  is  not  small,  lee  waves  generated  on 
the  ridge  convert  enough  energy  from  the  barotropic  tides  to  the  internal  tides  to  form 
tidal  bores  and  solitary  waves.  Increasing  stratification  at  ridge  depths  enhances  the 
generation  of  internal  waves,  particularly  at  the  diurnal  periods.  In  the  Luzon  Strait,  the 
slope  of  the  wave  beam  decreases  in  spring  and  summer  as  stratification  at  the  ndge 
depths  increases,  favoring  the  generation  of  internal  tides.  Without  the  presence  of  a 
strong  shallow  thermoeline,  internal  solitary  waves  are  not  observed  east  of  the  Luzon 
Strait.  In  the  northern  South  China  Sea,  internal  solitary  waves  are  likely  observed  from 
April  to  July  when  a  strong  shallow  thermoeline  is  present.  A  deep  mixed  layer  in  winter 
suppresses  the  production  of  internal  solitary  waves. 

Citation:  Shaw,  P.-T.,  D.  S.  Ko,  and  S.-Y.  Chao  (2009),  Inlemal  solitary  waves  induced  by  flow  over  a  ridge.  With  applications  io 
the  northern  Soulh  China  Sea,  J.  Geophys .  Res.,  114 ,  C02019,  doi:  10.1 029/2008JC005007. 


I.  Introduction 

[2]  Westward  propagating  internal  solitary  waves  (lSWs) 
in  the  northern  South  China  Sea  are  ubiquitous  in  satellite 
images  [c.g.,  Liu  et  ai ,  2004]  and  mooring  observations 
[Ramp  et  ai ,  2004].  In  the  satellite  images,  these  waves 
travel  toward  the  continental  margin  southeast  of  China  in 
the  form  of  long,  north  south  running  crests  [c.g.,  Liu  et 
ai ,  2004;  Zhao  et  ai ,  2004]  and  arc  reflected  and  dif¬ 
fracted  by  the  Tung-Sha  Island  [Chao  et  ai,  2006].  During 
the  passage  of  lSWs,  a  westward  velocity  of  2  m  s  1  and  a 
thermoeline  depression  of  150  m  have  been  recorded  in 
moored  observations  [Ramp  et  ai ,  2004;  Yang  et  ai ,  2004]. 
Most  ISWs  arc  first-mode  waves  with  a  depression  front 
but  second-mode  waves  have  been  observed  as  well  [Yang 
et  ai ,  2004],  The  propagation  speed,  estimated  to  be  1.5 
3.0  m  s  1 ,  is  close  to  the  phase  speed  of  the  first-mode 
internal  waves.  A  map  of  the  observed  ISWs  in  the 
northern  South  China  Sea  is  given  by  Zheng  et  al  [2007, 
Figure  2a] 

'Department  of  Marine,  Earth,  and  Atmospheric  Sciences,  North 
Carolina  Stale  Univcrsily,  Raleigh,  North  Carolina,  USA. 

’  Oceanography  Division,  Naval  Research  Laboratory,  Stcnnis  Space 
Center,  Mississippi,  USA. 

hlom  Poinl  Laboratory,  Center  for  Environmental  Science,  University 
of  Maryland,  Cambridge,  Maryland,  USA 

Copyright  2009  by  the  American  Geophysical  Union. 

0 1 48-0227/09/2008JC005007S09.00 


[3]  The  large  vertical  isopyenal  displacement  and  the 
strong  horizontal  How  induced  by  ISWs  are  of  concern  to 
martne  operations,  warranting  prediction  of  the  occurrence 
of  these  waves.  ISWs  at  many  locations  in  the  world 
ocean  are  generated  by  strong  barotropic  tidal  flow  over  a 
submarine  ridge  [c.g.,  Vlasenko  et  ai ,  2005].  In  the 
Luzon  Strait,  the  numerical  experiments  of  Niwa  and 
Hibiya  [2004]  show  that  M2  internal  tides  are  effectively 
generated  over  two  prominent  subsurface  ndges  by  the 
barotropic  tides  as  indicated  by  the  divergence  in  the 
integrated  internal  tidal  energy  flux.  Tracing  wave  rays  in 
satellite  images  clearly  indicates  that  the  Luzon  Strait  is 
the  source  region  of  these  waves  [Zheng  et  al .,  2007]. 
Using  phase  speed  estimated  from  climatological  stratifi¬ 
cation,  Zhao  and  Alford  [2006]  further  correlate  the 
arrival  of  ISW  packets  in  the  northern  South  China  Sea 
with  the  barotropic  tidal  flow  in  the  Luzon  Strait 
Amplitude  of  the  ISWs  has  also  been  found  to  correlate 
with  the  range  of  the  barotropic  tides  in  the  Luzon  Strait 
[Lien  et  ai ,  2005]. 

[4]  In  the  satellite  images,  ISWs  generated  in  the 
Luzon  Strait  propagate  westward  into  the  South  China 
Sea.  Few  or  no  ISWs  are  documented  over  the  ridges  and 
to  the  cast  of  the  Luzon  Strait  The  occurrence  of  ISWs 
also  depends  on  season  Statistical  analysis  shows  the 
ubiquitous  presence  of  ISWs  in  the  northern  South  China 
Sea  from  April  to  July  [Zheng  et  ai ,  2007].  If  the  same 
oceanic  conditions  arc  used  on  both  sides  of  the  L  uzon 
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Strait,  numerical  simulations  [e.g.,  Shaw  and  Chao ,  2006] 
indicate  that  ISWs  propagate  eastward  and  westward, 
unable  to  explain  the  lack  of  ISWs  east  of  the  Luzon 
Strait.  In  addition,  how  stratification  in  the  northern  South 
China  Sea  affects  the  seasonal  variation  of  ISWs  needs  to 
be  clarified 

[5]  Two  mechanisms  of  1SW  generation  over  a  ridge  have 
been  suggested  in  the  weakly  nonlinear  theory  [e.g., 
Helfrich  and  Melville ,  2006].  In  the  lee  wave  generation 
mechanism  [Maxworthy\  1979],  a  depression  in  the  thermo- 
elinc  is  first  produced  in  the  lee  of  a  ridge  by  the  barotropic 
tidal  flow.  As  the  tide  reverses,  the  depression  moves  away 
from  the  ridge  in  the  direction  of  the  flow  to  form  a  rank- 
ordered  solitary  wave  packet.  Apel  el  al.  [1985]  used  this 
mechanism  to  describe  the  ISW  generation  in  the  Sulu  Sea. 
The  second  mechanism  is  associated  with  internal  tide 
generation.  Fronts  develop  in  large-amplitude  internal  tides 
produced  by  flow  over  a  ndge.  Rank-ordered  wave  pack¬ 
ets  then  form  at  the  wavefronts  by  nonlinear  dispersion 
[Lee  and  Beardsley ,  1974].  Zhao  et  al.  [2004]  linked  the 
ISW  generation  in  the  northern  South  China  Sea  to  this 
mechanism  However,  observations  of  the  ISWs  in  the 
northern  South  China  Sea  are  unable  to  differentiate  the 
two  mechanisms. 

[6]  Most  ISW  theories  assume  weak  nonlinearity  and 
small  topography  [Vlasenko  et  al .,  2005;  Helfrich  and 
Melville ,  2006].  In  the  South  China  Sea,  the  thermocline 
depression  associated  with  the  ISWs  is  large,  and  ridges  in 
the  Luzon  Strait  occupy  a  good  fraction  of  the  water  depth. 
Beyond  weak  nonlinearity  and  small  topography,  this  study 
simulates  ridge-generated  internal  tides  and  their  transfor¬ 
mation  to  ISWs  in  a  nonhydrostatie  numerical  model  under 
idealized  conditions  This  work  is  part  of  a  larger  project 
sponsored  by  the  Office  of  Naval  Research.  The  project 
includes  field  observations,  remote  sensing,  theoretical 
studies,  and  numerical  efforts  with  a  goal  to  develop 
scientific  understating  that  leads  to  a  predictive  capability 
of  ISWs.  Toward  this  goal,  this  process  study  features  the 
oceanic  conditions  favoring  the  development  of  ISWs  in 
the  northern  South  China  Sea.  Restricting  our  goal  as  such, 
the  simulations  do  not  purposely  reproduce  the  quantitative 
aspects  of  ISW  properties,  e.g.,  the  observed  wave  ampli¬ 
tude  and  the  distance  from  the  ridge  where  nonlinear  ISWs 
first  appear. 

[7]  The  organization  of  this  paper  is  as  follows.  Section  2 
briefly  reviews  the  model  formulation.  Results  from  the 
simulation  are  given  in  Section  3.  From  the  simulation 
results,  optimal  conditions  for  the  generation  of  ISWs  are 
discussed  in  Section  4  and  summarized  in  Section  5.  It  is 
found  that  prediction  of  ISWs  requires  knowledge  of  the 
barotropic  tides  and  stratification  at  the  ridges  in  the  Luzon 
Strait  as  well  as  the  upper  ocean  stratification  in  the  northern 
South  China  Sea. 


equations  with  Boussinesq  and  rigid  lid  approximations. 
The  governing  equations  are 

—  4-  2Uk'  x  v  =  -  —  Vp  -  +  AV2,, v  -4  i'-—,  (1 ) 

Dt  p0  p0  <):- 


V  •  v  =  0,  and  (2) 


Dp 


(fp 


757  -  KXllP+k7^2  ' 


<?.-2 


(3) 


where  v  is  the  three-dimensional  velocity  vector  (//,  v,  w)  in 
the  (.v,  y,  z)  coordinate,  k'  is  a  unit  vector  pointing  upward 
from  the  North  Pole,  k  is  the  local  upward  unit  vector,  p  is  the 
density  after  removing  a  reference  value  po(=1028  kg  m  ■1)t 
fi  is  the  Earth's  rotation  rate,  and  p  is  the  pressure.  The 
eddy  viscosity  and  diffusivity  arc  constant,  so  that  A  K 
4  itT  s  *,//=!  x  10  4  m2  s  \  and  k  =  0. 1  x  10  4  m2  s  l. 
Equation  (1)  includes  both  the  horizontal  and  vertical 
components  of  the  Coriolis  acceleration,  and  the  traditional 
approximation  of  neglecting  the  vertical  Coriolis  accelera¬ 
tion  is  not  used. 

[9]  Tlie  finite  difference  form  of  the  governing  equations 
and  the  solution  procedure  are  described  by  Shaw  and  Chao 
[2006].  A  brief  outline  is  given  here.  The  Arakawa-C  grid 
system  is  used  for  the  spatial  discretization  of  equations  (1 ), 
(2),  and  (3).  Vertical  discretization  follows  the  z  coordinate 
formulation.  Spatial  derivatives  are  center  differenced  to 
second-order  accuracy.  Integration  in  time  follows  the 
Adams-Bashforth  scheme.  At  each  time  step,  the  baroclinie 
pressure,  barotropic  pressure,  and  nonhydrostatie  pressure 
are  computed  separately.  The  baroclinie  pressure  is  calcu¬ 
lated  by  vertically  integrating  the  nght-hand  side  of  the  z 
momentum  equation.  Setting  the  divergence  of  the  horizon¬ 
tal  transport  to  zero  then  produces  a  two-dimensional 
Poisson  equation  for  the  barotropie  pressure.  The  equation 
is  solved  by  matrix  inversion.  A  three-dimensional  Poisson 
equation  is  derived  for  the  nonhydrostatie  pressure  by  taking 
the  divergence  of  equation  (I)  and  setting  it  to  zero  using 
equation  (2).  The  equation  is  solved  by  iteration  using  the 
preconditioned  conjugate  gradient  method  [e.g.,  Pozrikidis , 
1997]. 

[10]  The  experiments  were  earned  out  in  a  channel  at 
20°N,  the  latitude  of  the  Luzon  Strait.  The  x  axis  is  toward 
east  along  the  channel  axis,  and  the  channel  depth  is  chosen 
to  be  800  m  for  simulation  of  internal  waves  in  the  upper 
ocean.  A  “bell-shaped”  ridge,  uniform  in  y  direction,  is 
placed  at  x  =  0  km.  The  ridge  height  h  is  given  by 


1  +  (x/L)2  ’ 


(4) 


2.  Model  Description 

[k]  The  three-dimensional  nonhydrostatie  numerical 
model  of  Shaw  and  Chao  [2006]  is  used  in  this  study.  Chao 
et  al.  [2006,  2007]  have  employed  this  model  in  earlier 
studies  of  ISWs  in  the  northern  South  China  Sea.  The  model 
uses  three-dimensional  momentum,  continuity,  and  density 


where  hm  and  L  are  the  maximum  height  and  the  scale  width 
of  the  ridge,  respectively.  In  this  study,  hm  =  600  m,  and  the 
top  of  the  ndge  is  200  m  below  sea  surface.  The  scale  width 
L  varies  from  1 5  km  to  60  km  (Table  I ).  The  west  and  east 
boundaries  are  located  at  x  =  -200  and  200  km, 
respectively.  Two  grid  cells  are  used  in  the  y  direction  A 
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Table  1.  List  of  Experiments  in  this  Study3 


Fxpcrimcnt 

L  (km) 

jV0  (rad  s  ') 

«o  (m  s  ') 

Phase  Speed  (ms  ') 

la 

15 

0015 

0.1 

1.52 

lb 

30 

0  015 

0.1 

1.52 

lc 

45 

0  015 

0.1 

1.52 

Id 

60 

0.015 

0.1 

1.52 

2a 

45 

0015 

0.2 

1.52 

2b 

60 

0.015 

0.2 

1.52 

3  a 

15 

0.008 

0.1 

0.76 

3b 

30 

0.008 

0.1 

0.76 

4a 

30 

0.0 15b 

0.1 

1.48 

4b 

30 

0  015" 

0.1 

1.05 

‘‘The  scale  width  of  the  ridge  (A)  is  defined  in  equation  (4),  and  the 
buoyancy  frequency  (jV0)  is  given  by  equation  (7).  The  amplitude  of 
barotropic  tides  is  m0.  The  last  column  is  the  phase  velocity  of  the  first  mode 
wave  in  deep  water  away  from  the  ridge. 

hMixcd  layer  (N  2  x  10  4  rad  s  ')  in  lhc  upper  50  m 
LMixcd  layer  (N  2  x  10  4  rad  s  ’)  in  the  upper  150  m. 


periodic  condition  in  y  results  in  solutions  with  no  y 
variation,  effectively  running  a  two-dimensional  simulation 
in  a  three-dimensional  model.  The  vertical  grid  size  is  10  m, 
and  the  horizontal  grid  size  is  200  m. 

[ i  t ]  The  density  field  consists  of  ambient  vertical  strati¬ 
fication  pa  (z)  and  the  perturbation  density  anomaly  p  (x,  z): 

/'  =  /'„(-)  +  ?'(*.-)•  (5) 

[12]  Except  for  the  mixed  layer,  the  ambient  stratification 
is  given  by 

=  ~~2  1  +  tanHC  /f)]’ 

where  A p  is  the  density  scale.  The  center  of  the  thermocline 
is  at  depth  z0(=— 50  m),  and  the  thickness  of  the  thermocline 
is  given  by  D(=  120  m).  The  buoyancy  frequency  is 

N  =  Hx/pM^'!2=  (0)  “Sh  1  (~) 

=  A/nCosh-1  (?) 

[13]  In  most  eases  (Experiments  1  and  2),  A p  is  chosen  to 
be  6  kg  m  3,  resulting  in  a  maximum  buoyancy  frequency 
N(z  =  z0)  =  iV0  =  0.015  rad  s_l.  Figure  la  shows  vertical 
profiles  of  density  and  buoyancy  frequency.  The  first  three 
normal  modes  for  horizontal  velocity  and  perturbation 
density,  calculated  by  solving  the  eigenvalue  problem 
numerically,  arc  plotted  in  Figures  lb  and  le,  respectively. 
The  zero  crossing  depth  for  the  mode  1  horizontal  velocity 
(also  the  depth  of  maximum  density  perturbation)  is  200  m. 
The  phase  speeds  of  the  first  three  baroelinie  modes  are 
1  52,  0.64  and  0.41  m  s  1  in  the  deep  water,  respectively. 
Over  the  ridge  top,  the  local  phase  speed  for  the  mode  1 
wave  is  calculated  to  be  0.89  m  s  1  using  the  local  water 
depth. 

[14]  Radiation  conditions  for  the  first-mode  baroelinie 
wave  are  adopted  at  the  open  boundaries  at  x  -  —200  and 
200  km.  A  rigid  lid  is  used  on  the  surface  for  computational 
efficiency.  Because  of  the  vast  difference  in  phase  speed 
(1.5  m  s  1  for  baroelinie  waves  versus  90  m  s  for 


barotropic  waves),  imposing  the  rigid  lid  approximation 
should  not  affect  the  amplitude  and  phase  speed  of  internal 
tides  over  the  ridge  significantly.  No  normal  flow  conditions 
are  used  on  all  solid  boundaries,  and  a  no-slip  condition 
dictates  the  tangential  velocity  on  the  bottom  Tidal  forcing 
is  included  by  specifying  the  sinusoidal  barotropic  inflow/ 
outflow'  velocity  on  the  two  x  boundaries.  The  tidal  current 
starts  from  zero  at  /  -  0  as  an  ebb  tide  with  eastward  flow.  In 
most  experiments,  the  maximum  tidal  velocity  on  open 
boundaries  (w0)  is  0  1  m  s  V  Although  the  semidiurnal 
tides  are  used  in  all  experiments  to  facilitate  discussion, 
similar  conclusions  can  be  extended  to  diurnal  tides.  All 
experiments  use  a  time  step  of  31.05  s,  which  is  1/120  of  the 
semidiurnal  tidal  period  (12.42  hours).  All  experiments  run 
for  five  tidal  cycles. 

[15]  Table  1  shows  four  sets  of  experiments  conducted  in 
this  study  The  first  set  (Experiments  la  to  Id  with  Exper¬ 
iment  lb  as  the  control  run)  examines  the  sensitivity  of  1SW 
generation  to  the  ridge  width  or  L.  The  second  set  (Experi¬ 
ments  2a  to  2b)  increases  the  barotropic  tidal  speed  twofold 
using  two  ridge  widths.  The  third  set  (Experiments  3a  to  3b) 
deviates  from  the  first  set  by  halving  the  buoyancy  frequency 
in  the  water  column  using  two  ridge  widths.  In  the  fourth  set 
(Experiments  4a  to  4b),  a  mixed  layer  is  imposed  at  the 
surface. 

3*  Results 

3.1.  Internal  Wave  Generation 

[16]  lSWs  arc  associated  with  sharp  depression  fronts, 
strong  downward  flow,  and  convergent  surface  currents. 
Among  them,  contour  plots  of  the  u  velocity  in  the 
numerical  result  give  the  best  visual  effects  to  identify 
lSWs.  Figure  2  shows  snapshots  of  u  velocity  contours 
between  2.5  and  3  tidal  cycles  after  initialization  for 
Experiment  lb,  the  control  run.  In  this  time  span,  the 
barotropic  tidal  current  is  to  the  left  or  negative.  In  the 
Luzon  Strait,  this  would  be  the  flood  phase  of  the  semidi¬ 
urnal  tidal  current.  As  the  flood  tide  intensifying  from  31.05 

f>B  (kg/m3)  p  p 


Figure  1.  (a)  Density  pa  in  kg  m  and  buoyancy 
frequency  N  in  rad  s  1  in  Experiments  1  and  2.  The 
corresponding  normal  modes  are  shown  (b)  for  horizontal 
velocity  and  (e)  for  density  perturbation.  The  phase  speeds 
are  1.52,  0.64,  and  0.41  ms  1  for  the  first  three  modes. 
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Figure  2.  Contour  plots  of  u  velocity  in  m  s  1  for  Experiment  lb  at  (a)  31.05,  (b)  33. 12,  (e)  35.19,  and 
(d)  37.26  hours  during  the  flood  tide  (barotropie  tidal  flow  to  the  left)  The  contour  interval  is  0.05  ms 


hours  to  33.12  hours,  two  slanting  wave  beams  of  strong 
westward  flow  extend  upward  from  a  point  over  the  ridge 
top.  In  the  bands,  eontour  lines,  whieh  arc  also  constant 
phase  lines,  move  downward  from  33.12  hours  to  35.19 
hours,  indieating  upward  energy  propagation  along  internal 
wave  beams.  Some  energy  also  propagates  downward  from 
the  ridge  top  along  both  faees  of  the  ridge.  Similar  wave 
beams  are  displayed  in  the  internal  wavefield  generated  by 
an  oscillating  point  souree  in  the  classical  laboratory  experi¬ 
ments  of  Mowbray  and  Rarity  [1967;  see  Gill ,  1982,  p. 
137].  Thus,  internal  waves  in  Figure  2  are  produced  by  the 
tidal  current  over  the  top  of  the  ridge.  Reflection  of  the 
wave  beams  by  the  surface  at  x  =  —  1 5  km  is  clearly  seen  at 
35.19  hours  (Figure  2e).  The  stronger  one  of  the  two 
upward  wave  beams  is  on  the  downstream  side  of  the  ridge. 
The  partiele  speeds  on  top  of  the  ridge  in  Figure  2  are  less 
than  0.5  ms1,  which  gives  a  moderate  local  Froude 
number  about  0.5  for  a  phase  speed  of  0.89  ms1.  Vlasenko 
et  al.  [2005]  concluded  that  nonlinear  e  fleets  do  not  play  a 
dominant  role  in  the  generation  of  internal  tides  over 
inclined  bottom  topography  at  low  Froude  numbers. 

[ 1 7]  The  two  wave  beams,  reflected  downward  by  the  sea 
surface,  are  refracted  upward  at  x  -  —30  km  as  they  enter 
deep  water  with  decreasing  stratification  After  zigzagging, 
the  upper  ocean  waveguide  transforms  the  wave  to  hori¬ 
zontally  propagating  shallow  water  waves  in  the  hydrostatic 
regime.  For  example,  the  first-mode  baroelinie  wave  struc¬ 
ture  is  established  on  the  left  side  of  the  ridge  after  the 
sceond  reflection  from  the  surface  (to  the  left  of  x  =  -60  km 
in  Figure  2a).  Intensification  of  wavefronts  occurs  quickly 


thereafter.  For  example,  a  tidal  bore  is  clearly  seen  at  jc 
—  71  km  in  Figure  2a.  Six  hours  later,  the  tidal  bore 
develops  into  a  rank-ordered,  internal  solitary  wave  train 
at  jc  =  -110  km  (Figure  2d).  To  its  left  is  another  tidal  bore 
generated  one  tidal  cycle  earlier. 

[is]  On  the  left  side  of  the  ridge,  a  tidal  bore  forms  on  the 
leading  edge  of  a  negative  velocity  region.  At  generation, 
negative  u  velocity  is  associated  with  the  flood  tides.  Thus, 
the  two  wave  trains  on  the  left  side  of  the  ridge  in  Figure  2 
were  generated  when  the  barotropie  flow  on  the  ridge  top 
ehanged  from  ebb  to  flood  tides  during  the  first  two  tidal 
eyeles.  On  the  right  side  of  the  ndge,  the  corresponding  tidal 
bores  appear  on  the  leading  edge  of  a  region  of  positive  u 
velocity  and  are  associated  with  the  beginning  of  ebb  tides 
at  generation. 

[19]  Contour  plots  of  perturbation  density  provide  more 
insights  into  the  wave  generation  mechanism  (Figure  3).  At 
the  slack  water  before  the  flood  tide  (hour  31  05),  the 
isopyenal  has  been  depressed  downward  on  the  right  (lee) 
side  of  the  ridge  by  an  earlier  ebb  tide  (Figure  3a).  Above 
the  ndge,  isopycnals  are  depressed  along  a  path  toward  left 
following  a  wave  beam  in  Figure  2a,  As  the  barotropie  tide 
turns  toward  left  (Figure  3b),  the  signal  in  isopyenal 
depression  is  refleeted  by  the  surfaee  as  in  Figure  2.  The 
first-mode  structure  in  isopyenal  depression  appears  at 
37.26  hours  (Figure  3d).  This  depression  will  eventually 
intensify  similar  to  the  one  at  x  =  -94  km  produced  one 
tidal  cycle  earlier  in  Figure  3e.  Figure  3  shows  that 
westward  propagating  internal  tidal  bores  originate  from 
lee  waves  produced  by  the  ebb  tide  on  the  east  side  of  the 
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Figure  3.  Contour  plot  of  perturbation  density  (p)  in  kg  m  for  Experiment  1  b  at  (a)  3 1 .05,  (b)  33. 1 2, 
(c)  35.19,  and  (d)  37.26  hours  during  the  flood  tide  The  contour  interval  is  0.5  kg  m  3  The  maximum 
velocity  is  0.5  m  s  1  for  u  and  0.01  ms  1  for  w. 


ridge.  Tidal  bore  generation  can  be  traced  back  to  the  start 
of  a  flood  tide,  in  agreement  with  the  timing  inferred  from 
Figure  2. 

3.2.  Formation  of  Internal  Solitary  Waves 
[20]  In  Experiment  lb,  lSWs  eventually  develop  in  a  tidal 
bore,  c.g.,  at  jc  =  —180  km  and  110  km  in  Figure  2d. 
Contours  of  perturbation  density  p  and  the  velocity  vectors 
associated  with  a  well-developed  solitary  wave  train  at  hour 
43.37  from  Experiment  lb  are  shown  in  Figure  4.  The 
leading  depression  is  at  jc  =  —142  km  with  an  amplitude  of 
30  m,  followed  by  a  smaller  depression  at  jc  =  —138  km. 
The  tidal  phase  corresponds  to  the  slack  water  before  the 
flood  tide  The  barotropic  flow  vanishes  at  this  time,  and  the 
(w,  w)  velocity  vector  represents  the  flow  field  of  the  solitary 
wave.  The  structure  is  described  by  the  first  baroclinie  mode 
(Figure  1).  The  maximum  speed  at  the  surface  is  0,4  ms  1 
and  is  located  at  x  142  km,  the  position  of  the  maximum 
depression.  The  distance  between  the  first  two  wave  troughs 
is  3.6  km.  Strong  convergence  produces  downward  motion 
about  0.6  km  ahead  of  the  maximum  horizontal  velocity  at 


the  front,  followed  by  alternating  upward  and  downward 
motion.  The  maximum  downward  velocity  is  0.01  ms  1 
The  magnitude  decreases  rapidly  behind  the  front. 

[21]  Propagation  of  internal  waves  at  z  —  —  5  m  is  shown 
in  the  contour  plot  of  the  u  velocity  on  the  x-t  plane 
(Figure  5).  Patches  of  minima  and  maxima  in  the  u  velocity 
contours  along  parallel  slanting  lines  show  outward  wave 
propagation  from  the  ridge  top  at  jc  =  0  km.  The  horizontal 
bands  arc  produced  by  the  alternating  barotropic  tidal 
current  Since  the  data  are  stored  every  hour,  the  crowded 
contour  lines  at  a  wavefront  cannot  be  fully  resolved, 
producing  small,  wavy  patches.  On  the  left  (right)  side  of 
the  ridge,  the  leading  edge  of  a  region  of  negative  (positive) 
u  velocity  intensifies  to  form  a  bore.  Propagation  of  lSWs 
on  the  left  (right)  side  of  the  ridge  can  be  traced  back  to  the 
beginning  of  the  flood  (ebb)  phase  of  the  barotropic  tides 
over  the  ridge  top.  The  timing  is  consistent  with  the  analysis 
in  Section  3.1  and  that  by  Zhao  and  Alford  [2006].  Esti¬ 
mated  from  the  slope  of  the  position  of  each  wavefront  on 
the  jc-/  plane,  the  propagation  speed  of  the  internal  bores  is 
fairly  constant  at  1 .65  ±  0.07  m  s  1 ,  or  about  6  km  h  1 .  This 
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Figure  4.  Perturbation  density  contours  and  the  velocity  vector  in  the  x-z  plane,  showing  a  well-developed 
solitary  wave  train  at  43.47  hours  in  Experiment  lb.  Contour  values  arc  from  —4  to  —0.5  kg  m  3  at  a 
contour  interval  of  0.5  kg  m  3  The  maximum  velocity  is  0.4  m  s  1  for  u  and  0.01  m  s  1  for  u\ 


value  is  higher  than  the  theoretical  value  of  1.52  m  s  !  for 
the  first  baroclinic  waves.  Nonlinear  effects  clearly  increase 
the  speed  of  wave  propagation. 

3.3.  Effect  of  the  Ridge  Width 

[22]  Figure  6  shows  contour  plots  of  the  m- velocity  at  the 

end  of  the  third  tidal  cycle  (/  =  37.26  hours)  for  ridges  with 
dilTcrcnt  width  scales  (Experiments  la,  lb,  Ic,  and  Id). 
Internal  wave  beams  exist  near  the  ridge  in  all  panels,  but 
the  wave  strength  decreases  with  increasing  ridge  width.  In 
all  cases,  vertically  propagating  waves  evolve  to  first-mode 
shallow  water  waves  about  100  km  away  from  the  ridge 
because  of  wave  reflection  and  refraction.  The  first-mode 
waves  then  intensify  to  form  ISWs  in  Figures  6a  and  6b. 
When  /  increases  to  45  km  (Figure  6c),  internal  waves  still 
emit  from  the  ridge  but  the  wavefront  at  x  =  —100  km 
shows  only  slight  intensification.  At  x  =  170  km,  the 

trough  produced  one  tidal  cycle  earlier  is  strongly  dissipated 
(Figure  6c).  Little  internal  wave  energy  leaves  the  ridge 
when  L  is  further  increased  to  60  km  (Figure  6d).  Thus,  the 
strength  of  the  baroclinic  tides  depends  on  the  ridge 
steepness.  Energy  is  converted  from  the  barotropic  tide  to 
the  baroclinic  tide  if  the  ridge  steepness  exceeds  a  critical 
value. 

[23]  Table  2  shows  the  first  three  local  extremes  in  u 
velocity  associated  with  left-propagating  waves  in  Figure  6. 
The  peak  velocities  at  the  two  tidal  bores  (x  -  —180  and 

1 10  km)  arc  given  by  u  \  and  w3,  while  u2  is  the  maximum 
positive  velocity  in  between  (jc  =  —130  km)  These  values 
describe  the  amplitude  of  the  baroclinic  wave  at  t  =  37.26 
hours  when  the  barotropic  tidal  velocity  is  zero.  In  partic¬ 
ular,  w,  and  m3  arc  the  troughs  in  the  baroclinic  How 
generated  during  the  first  two  flood  phases  of  the  barotropic 
tides,  and  u2  is  the  peak  speed  generated  during  the  second 
ebb  tide  (Figure  5).  The  last  column  indicates  the  presence 
of  ISWs  in  the  experiment.  When  ISWs  form  (Experiments 


la  and  lb),  the  value  of  u2  is  higher  than  the  maximum 
barotropic  flow  (0. 1  m  s  *).  Furthermore,  m,  and  m3  are 
much  greater  than  u2  because  of  intensification  of  the 
wavefronts.  Experiment  lc  is  near  the  borderline  of  ISW 
formation.  With  a  weak  front  at  jc  =  — 100  km,  u2  and  1/3  are 
about  the  same  as  the  barotropic  tidal  velocity,  but  M|  is 
smaller.  In  Experiment  Id,  the  baroclinic  velocity  u2  is  0.06 
ms  1  or  2/3  of  the  barotropic  velocity,  while  U\  and  m3  arc 
about  half  of  u2.  In  this  case,  no  ISWs  form.  Multiple 
troughs  arc  usually  associated  with  the  ISWs,  e  g.,  between 
x  =  180  and  —190  km  in  Figure  6a.  The  average  wave¬ 

lengths,  estimated  from  the  distance  of  the  first  two  com 
secutivc  troughs,  arc  4  km  and  3  km,  for  u  j  and  M3, 
respectively.  Similar  values  have  been  reported  in  the 
observed  ISWs  [c.g..  Ramp  ct  al.y  2004]  Thus,  surface 


u  (m/s)  at  z  =  -5  m 
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Figure  5.  Contour  plot  of  u  velocity  at  z  =  —  5  m  on  the  x-l 
plane  for  Experiment  lb  The  contour  interval  is  0.1  ni  s  1 . 
The  ridge  peak  is  at  jc  =  0  km 
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Figure  6.  Contour  plots  of  u  velocity  at  slack  water  before  Hood  tide  at  t  37.26  hours  for  four  scale 
widths  of  the  ridge  ( L ):  (a)  1 5  km,  (b)  30  km,  (c)  45  km,  and  (d)  60  km.  The  contour  interval  is  0.05  ms  1 . 


convergence  with  velocity  much  greater  than  the  maximum 
barotropie  tidal  current  is  a  telling  sign  of  the  presence  of 
ISWs. 

[24]  If  the  barotropie  velocity  in  Experiment  1  e  is  doubled 
to  0.2  ms1,  sharp  fronts  form  in  Experiment  2a  (Figure  7a). 
Compared  to  the  values  in  Experiment  Ie,  u2  is  nearly 
doubled  (0.19  ms  l),  but  U\  and  w3  are  three  times  that  in 
Experiment  1c  (Table  2),  indicating  1SW  formation.  When 
the  barotropie  velocity  in  Experiment  Id  is  doubled,  wt,  m2, 
and  w-<  in  Experiment  2b  increase  proportionally  with  all 
three  values  much  lower  than  the  barotropie  How  (Table  2). 
No  wavefronts  form  in  this  case  (Figure  7b).  In  Experiments 
2a  and  2b,  the  maximum  barotropie  tidal  velocity  is  0.8  ms  1 
on  top  of  the  ridge.  With  a  local  phase  speed  of  0.89  m  s  \ 
the  local  Froude  number  is  close  to  I .  Increase  in  the  strength 
of  lee  waves  as  the  Froude  number  approaching  I  may 
explain  why  ISWs  appear  in  Experiment  2a,  in  which  the 
ridge  slope  is  near  the  critical  value.  However,  for  a  slope 
below  the  critical  value  as  in  Experiment  2b,  ISWs  are  not 
generated  even  the  Froude  number  is  near  I . 

3.4.  Stratification 

[25]  In  the  linear  internal  wave  theory,  the  slope  of  the 
wave  beam  in  an  ocean  with  constant  buoyancy  frequency 
No  is 

tancf>  =  yj(u.2  ~P)/(Nl  -  u^)  (8) 

where  <f)  is  the  angle  between  the  wave  beam  and  the 
horizontal  axis,  tu  is  the  wave  frequency,  and  /  is  the 


Conolis  parameter  [e.g..  Gill,  1982,  p.  259],  Thus,  the 
critical  ridge  slope  is  inversely  proportional  to  N0  at  tidal 
frequencies.  Experiments  3a  and  3b  parallel  Experiments  la 
and  lb,  respectively,  except  for  a  twofold  reduction  of  the 
buoyancy  frequency  throughout  the  water  column  From 
equation  (8),  the  critical  slope  in  Experiments  3a  and  3b  is 
twice  that  in  Experiments  I  and  2.  The  phase  speeds  for  the 
first  two  normal  modes  are  reduced  to  0  76  and  0.32  m  s  \ 
respectively,  and  the  Froude  number  is  about  I  for  mode  1 
waves  in  both  Experiment  3a  and  3b.  The  contour  plot  of  u 
velocity  in  Experiment  3a  shows  multiple  troughs  trailing 
the  wavefront  (Figure  7e).  In  Table  2,  the  speeds  at  the  crest 
(«!  and  w3)  are  more  than  twice  that  at  the  trough  (m2), 
which  is  slightly  higher  than  the  barotropie  tidal  speed  (m0). 


Table  2.  First  Three  Local  Extremes  in  the  Velocity  of  the  Left* 
Propagating  Waves  at  t  =  37.26  ha 


Experiment  u0 

(m  s 

')  "1  (m  s  ') 

i/2  (ms1) 

u\  (ms1) 

ISW? 

la 

0.1 

-036 

0.15 

-0.40 

Yes 

lb 

0.1 

-0  29 

0.15 

-0.31 

Yes 

Ic 

0.1 

-0.08 

0.11 

-0.12 

No 

Id 

0.1 

-0.04 

006 

-0.03 

No 

2a 

0.2 

0.29 

0.19 

0.37 

Yes 

2b 

0.2 

-0.10 

0.13 

0.07 

No 

3a 

0.1 

-0.27 

0.13 

0.32 

Yes 

3b 

0.1 

-005 

0.06 

-0.03 

No 

4a 

0.1 

-0.22 

0.14 

-0.24 

Yes 

4b 

0.1 

Unidentifiable  Unidentifiable  Unidentifiable 

No 

“The  maximum  barotropie  tidal  flow  is  The  baroclinie  velocities  arc 
given  by  M|,  w2,  and  w*.  The  last  column  indicates  whether  1SW  is  generated 
in  each  experiment. 
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Figure  7.  Contour  plots  of  u  velocity  at  slack  water  before  Hood  tide  at  t  =  37.26  hours  for 
(a)  Experiment  2a,  (b)  Experiment  2b,  (e)  Experiment  3a,  (d)  Experiment  3b,  (e)  Experiment  4a,  and 
(0  Experiment  4b.  The  contour  interval  is  0.05  ms  1  in  all  eases  except  0.02  m  s  1  in  Figures  7d  and  7f. 


ISWs  are  generated  in  this  ease.  The  wave  amplitude  is  only 
slightly  lower  than  that  in  Experiment  la  (Table  2). 

[26]  If  L  in  Experiment  3a  is  doubled  (Experiment  3b), 
only  weak  internal  waves  propagate  away  from  the  ridge 
(Figure  7d).  The  internal  waves  are  so  weak  that  ISWs  no 
longer  develop.  The  u  velocity  in  Experiment  3b  is  quali¬ 
tatively  different  from  that  in  experiment  lb  The  speed  at 
the  trough  in  Experiment  3b  is  lower  than  the  barotropie 
tidal  speed,  indicating  no  wavefront  formation  (Table  2). 
The  speed  at  the  trough  is  also  lower  than  that  at  the  crest  as 
expected.  Again,  the  ndge  slope  must  exeeed  a  threshold  to 
convert  sufficient  energy  from  the  barotropie  tides  to 
internal  tides  for  1SW  formation  even  the  Froude  number 
is  close  to  1 . 

[27]  Experiments  4a  and  4b  use  the  same  stratification  as 
in  Experiment  lb  but  the  upper  ocean  now  contains  a  mixed 
layer  50  and  150  m  deep,  respectively.  The  mixed  layer 
is  weakly  stratified  with  a  buoyancy  frequency  of  2  x 
10  4  rad  s  \  allowing  internal  waves  at  the  semidiurnal 
frequency  to  exist  in  the  mixed  layer.  The  first-mode  phase 
speeds  decrease  to  1.48  m  s  1  and  1.05  m  s  1  in  Experi¬ 
ments  4a  and  4b,  respectively.  In  Figure  7e,  solitary  wave 
trains  appear  in  the  u  velocity  field  at  hour  37.26;  develop¬ 
ment  of  ISWs  is  not  affected  by  the  presence  of  a  shallow 
mixed  layer.  The  extreme  velocities  in  Table  2  show  some 
reduction  in  peak  velocities  ux  and  1/3  though.  If  the  mixed 
layer  depth  increases  to  1 50  m,  the  mixed  layer  effectively 
erases  the  thermoeline  at  z  =  -50  m  shown  in  Figure  la. 
Internal  wave  propagation  is  hardly  seen  (Figure  70-  Another 


experiment  has  been  conducted  to  further  distinguish  the 
role  of  the  thermoeline  in  ISW  generation.  The  stratification 
in  the  top  300  m  of  the  water  column  in  Experiment  lb  is 
replaced  by  one  with  constant  buoyancy  frequency  (N  = 
0  015  rad  s  ')  The  thermoeline  is  diffused  over  the  top 
300  m,  but  the  maximum  buoyancy  frequency  is  unchanged. 
The  result  (not  shown)  contains  only  linear  hydrostatic 
internal  waves  propagating  away  from  the  ndge,  pointing 
out  the  need  of  a  sharp  thermoeline  to  trap  wave  beams  for 
ISW  generation 

4.  Discussion 

[2«]  Barotropie  tides  over  a  ndge  generate  internal  tides 
propagating  away  from  the  ridge  top  in  slanting  wave 
beams.  Away  from  the  ndge,  an  upper  ocean  waveguide 
at  the  buoyancy  frequency  maximum  (Figure  la)  traps 
internal  waves  because  of  wave  reflection  and  refraction 
Wave  propagation  eventually  becomes  horizontal  some 
distance  away  from  the  ridge.  The  strength  of  the  internal 
tides  is  crucial  to  the  development  of  tidal  bores  and  ISWs 
Significant  energy  is  converted  from  the  barotropie  tides  to 
the  internal  tides  only  when  the  ridge  slope  exceeds  a 
critical  value.  Thus,  a  steep  ndge  would  be  most  efficient 
in  generating  internal  tides  in  the  ocean.  Linear  theory  [e.g.. 
Hurley  and  Keady\  1997]  and  laboratory  experiments  [e.g., 
Zhang  et  al ,  2007]  suggest  that  internal  wave  generation  is 
possible  when  the  ridge  slope  exceeds  the  slope  of  the  wave 
beam.  The  numerical  expenments  in  section  3  validate  this 
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criterion  when  the  flow  is  nonlinear  and  the  ridge  height  is  a 
large  fraction  of  the  water  depth.  If  the  ridge  slope  is  less 
than  the  critical  value,  only  weak  internal  waves  are 
produced. 

[29]  Tidal  bores,  if  formed,  are  always  located  on  the 
leading  edge  of  the  down-wave  surface  flow.  In  the  northern 
South  China  Sea,  westward  propagating  ISWs  are  associated 
with  the  leading  edge  of  the  westward  flow  in  the  upper 
ocean  At  generation,  the  westward  surface  baroclinic  flow 
occurs  when  the  barotropic  tidal  flow  over  the  ndge  changes 
from  ebb  (eastward)  to  flood  (westward).  The  analysis  by 
Zhao  and  Alford  [2006]  shows  that  ISWs  in  the  northern 
South  China  Sea  are  associated  with  westward  barotropic 
tides  in  the  Luzon  Strait.  The  numerical  result  is  consistent 
with  the  observed  wave  behavior. 

[30]  Equation  (8)  gives  the  critical  slope  as  a  function  of 
wave  frequency  and  buoyancy  frequency.  The  wave  beam 
becomes  more  horizontal  as  the  buoyancy  frequency 
increases.  Thus,  strong  stratification  at  ridge  depths  favors 
generation  of  baroclinic  tides.  Furthermore,  the  internal 
wave  beam  is  more  horizontal  at  the  diurnal  period  than 
at  the  semidiurnal  period  under  the  same  stratification  A 
ndge  is  therefore  more  likely  to  become  supercritical  for 
diurnal  tides  than  for  semidiurnal  tides  although  simulations 
of  diurnal  tides  are  not  explicitly  illustrated  in  this  study. 
Ramp  et  al.  [2004]  observed  strong  ISWs  at  the  diurnal 
period  in  this  region  (their  24-hour  or  type  a  waves). 
Diurnal  barotropic  tides  are  stronger  than  the  semidiurnal 
tides  in  the  L  uzon  Strait,  but  a  supercritical  ndge  slope  for 
the  diurnal  tides  may  also  contribute  to  the  difference  in 
amplitude  of  internal  waves. 

[31]  The  discussion  based  on  linear  theory  holds  when  the 
local  Froude  number  on  the  ridge  top  is  much  less  than  I  In 
the  linear  wave  limit,  two  symmetrical  wave  beams  emit 
from  the  ndge  top.  Lee  waves  arc  produced  over  a  large- 
amplitude  obstacle  by  barotropic  oscillatory  flow  in  the 
Froude  number  range  between  0.5  and  1  75  in  the  labora¬ 
tory  and  numerical  expenments  of  Matsuura  and  Hihiya 
[1900].  Westward  propagation  of  depressed  isopycnals  in 
Figure  3  (Froude  number  ~  0.5)  is  consistent  with  the  lee 
wave  theory.  During  the  ebb  tide,  isopycnals  are  depressed 
on  the  east  side  of  the  ridge  because  of  energy  accumulation 
by  internal  waves  propagating  against  the  eastward  tidal 
current.  As  the  tidal  current  reverses  to  the  west,  the 
depression  propagates  westward  over  the  ndge  to  enhance 
the  wave  beam  on  the  left  side  of  the  ridge.  Large  amplitude 
internal  waves  occur  at  Froude  number  ~  l  over  a  narrow 
ridge  in  the  experiments  of  Matsuura  and  Hihiya  [1990]. 
The  Froude  number  is  affected  by  the  strength  of  the  tidal 
current  and  the  local  phase  velocity  over  the  ridge  top.  Zhao 
and  Alford  [2006]  rejected  the  lee  wave  mechanism  on  the 
basis  of  the  lack  of  eastward  flow  immediately  before  the 
westward  flow  that  produces  ISWs.  However,  by  examining 
their  Figure  8,  a  period  of  deceleration  and  acceleration 
always  appears  in  the  westward  current  even  if  two  ISW- 
producing  peaks  in  the  westward  current  are  not  separated 
by  eastward  flow.  It  is  likely  that  two  consecutive  westward 
bursts  of  lee  waves  arc  released  in  one  diurnal  eyele  as  a 
response  to  mixed  semidiurnal  and  diurnal  tidal  forcing.  The 
lee  wave  bursts,  originating  from  a  depression  produced  on 
the  cast  side  of  the  ridge  by  the  previous  eastward  tidal  flow, 


arc  separated  by  decelerating  and  accelerating  tides  rather 
than  eastward  flow. 

[32]  It  seems  that  both  lee  waves  and  internal  tides  play  a 
role  in  the  ISW  generation.  The  ridge  slope  controls  the 
generation  of  linear  internal  tides  from  a  ridge  but  energy 
needed  to  form  nonlinear  ISWs  has  to  come  from  lee  waves. 
Thus,  ISW  generation  requires  a  supercritical  slope  as  well 
as  a  Froude  number  close  to  1  over  the  ridge  top.  For  a  ridge 
slope  greater  than  the  critical  value,  internal  waves  are 
generated  and  may  evolve  into  nonlinear  ISWs  if  the  Froude 
number  is  not  too  small.  The  local  Froude  number  does  not 
constitute  an  overbearing  criterion  for  ISW  generation 
though.  ISWs  are  not  generated  when  the  ridge  slope  is 
less  than  the  critical  value  even  when  the  local  Froude 
number  is  close  to  one. 

[33]  Development  of  ISWs  involves  not  only  the  conver¬ 
sion  of  energy  from  the  barotropic  tides  to  the  baroclinic 
tides  but  also  the  formation  of  horizontally  propagating 
internal  waves  at  the  thermoclinc  [ Gerkema ,  2001 ;  Akylas  el 
al.,  2007]  Thus,  generation  of  ISWs  is  a  two-step  process. 
First,  internal  wave  beams  arc  produced  by  the  barotropic 
tides  over  a  ndge.  This  step  requires  a  supercritical  ridge 
slope  and  a  Froude  number  that  is  not  too  small.  In  the 
second  step,  an  upper  ocean  waveguide  produced  by  a  sharp 
thermoclinc  is  needed  to  trap  internal  wave  energy.  L.arge- 
amplitudc  internal  tides  are  produced  if  the  thcmiocline  is 
sharp;  subsequent  nonlinear  intensification  in  the  internal 
tides  leads  to  the  formation  of  tidal  bores  and  ISWs. 
Without  a  strong  thermoclinc,  only  linear  internal  waves 
propagate  away  from  the  ridge  A  shallow  upper  ocean 
mixed  layer  does  not  affect  the  formation  of  ISWs,  but  a 
deep  mixed  layer  that  erases  the  sharp  thermoclinc  would 
prevent  the  development  of  ISWs. 

[34]  The  Kuroshio,  flowing  northward  along  a  path  east 
of  the  Luzon  Strait,  is  a  boundary  between  the  warm  water 
in  the  Philippine  Sea  and  the  colder  subsurface  water  in  the 
South  China  Sea.  The  thermoclinc  is  noticeably  deeper  on 
the  cast  side  of  the  Kuroshio  than  in  the  South  China  Sea. 
The  position  of  the  Kuroshio  affects  the  strati lleation  at 
ndge  depths  in  the  Luzon  Strait  and  thus  internal  wave 
generation.  The  surface  position  of  the  Kuroshio  has  been 
studied  [e.g..  Yuan  et  al.,  2006],  but  few  observations  on  the 
variation  of  its  subsurface  boundary  arc  available.  There¬ 
fore,  the  climatology  data  arc  used  to  describe  the  change  in 
stratification.  The  monthly  buoyancy  frequency  distribution 
at  400  m  depth  near  the  Luzon  Strait  in  Figure  8  is 
constructed  from  the  Navy  Ocean  Climatology  Database 
(GDEM  version  3.0),  which  is  compiled  from  observed 
temperature  and  salinity  profiles.  The  buoyancy  frequency 
increases  from  the  South  China  Sea  to  the  Philippine  Sea, 
ranging  from  0.0061  rad  s  1  (3.5  cycles  h  l)  in  dark  blue  to 
0.0079  rad  s  1  (4.5  cycles  h  l)  in  dark  red.  Two  ridges  are 
present  in  the  Luzon  Strait,  the  east  ndge  following  the 
island  chain  to  the  east  and  the  west  ridge  extending 
southward  from  the  southern  tip  of  Taiwan.  In  Figure  8, 
stratification  at  ridge  depths  is  stronger  from  March  to 
August  than  the  rest  of  the  year.  According  to  the  numerical 
simulation  in  Section  3,  spring  and  summer  should  be  the 
favorite  period  for  the  generation  of  internal  wave  beams. 

[35]  Figure  9,  also  from  the  Navy  Ocean  Climatology 
Database,  shows  mean  vertical  profiles  of  density  and 
buoyancy  frequency  on  the  two  sides  of  the  Luzon  Strait 
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Figure  8.  Monthly  mean  buoyancy  frequency  at  400  m  in  the  Lu/on  Strait.  The  color  varies  progressively 
from  dark  blue  at  3.5  eyeles  h  (0.0061  rad  s  ')  to  dark  red  at  4.5  eyeles  h  1  (0.0079  rad  s  ').  The  islands 
at  top  and  bottom  of  Figure  8  are  Taiwan  and  Luzon,  respectively. 


between  Apnl  and  June.  The  left  panel  indicates  that  the 
thermoeline  is  stronger  and  shallower  on  the  west  side  than 
on  the  east  side  of  the  Luzon  Strait.  The  buoyancy  frequen¬ 
cy  profile  in  the  right  panel  shows  higher  values  below  150 
m  on  the  east  side  of  the  Luzon  Strait  than  in  the  South 
China  Sea,  consistent  with  Figure  8.  Because  the  steeper, 
taller  east  ridge  is  more  likely  in  waters  of  stronger 
stratification  than  the  west  ridge;  the  east  ridge  is  likely 
the  major  source  of  ISWs  in  the  northern  South  China  Sea. 
The  buoyancy  frequency  profiles  also  show  a  shallow  layer 
of  high  buoyancy  frequency  in  the  South  China  Sea  but  not 
east  of  the  Luzon  Strait.  Missing  an  upper  ocean  waveguide 
associated  with  a  strong,  shallow  thermoeline  may  explain 
the  lack  of  ISWs  east  of  the  Luzon  Strait. 

[ja]  The  shallow  thermoeline  in  the  northern  South  China 
Sea  favors  ISW  formation,  a  deep  mixed  layer  weakens  the 
thermoeline  and  thus  hinders  the  ISW  formation.  Local 
cooling  and  wind-induced  mixing  may  produce  a  mixed 
layer,  but  intrusion  of  surface  water  from  the  Philippine  Sea 
in  winter  [Shaw,  1991;  Yuan  et  al .,  2006]  effectively  erases 
the  sharp  thermoeline.  Qu  et  al.  [2007,  Figure  10]  calculated 
the  monthly  depth  of  the  isothermal  surfaee  that  coincides 
with  the  base  of  the  mixed  layer  in  the  South  China  Sea.  In  the 


region  west  of  Luzon  Strait,  the  mixed  layer  depth  decreases 
from  January  to  June  and  then  increases  to  a  maximum  in 
November.  The  shallowest  mixed  layer  occurs  from  April  to 


Figure  9.  Vertical  profiles  of  mean  density  (left)  and 
buoyancy  frequency  (right)  averaged  over  April  June  at 
20°N,  1 1 8°E  (solid  line)  and  at  20°N,  124°E  (dot  line). 
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August,  a  favorite  period  for  ISW  formation  as  predicted  by 
the  present  study  and  in  the  study  of  Zhao  and  Alford  [2006]. 

5.  Conclusions 

[37]  The  numerical  experiments  in  this  study  identify  the 
process  crucial  to  the  formation  of  ISWs  in  the  northern 
South  China  Sea.  The  model  can  be  fine-tuned  to  match  the 
observed  parameters,  e.g.,  the  amplitude  of  ISWs,  the 
distance  where  nonlinear  ISWs  begin  to  form,  etc.,  but  it 
is  not  the  intention  of  this  paper.  Rather,  ISW  generation  by 
the  barotropie  tides  is  demonstrated  in  the  simplest  model 
configuration.  The  oceanographic  conditions  that  may  affect 
ISW  generation  arc  identified  and  arc  applied  to  the 
northern  South  China  Sea. 

[ix]  Generation  of  ISWs  is  a  two-step  process.  First, 
internal  wave  beams  arc  produced  in  the  Luzon  Strait  by 
the  barotropie  ttdes  over  ridges  with  supercritical  slopes. 
The  critical  value  for  slope  is  determined  by  stratification, 
which  is  in  turn  affected  by  the  position  of  the  Kuroshio. 
The  stronger  stratification  at  400  m  in  spring  and  early 
summer  near  the  ridges  in  the  Luzon  Strait  may  enhance  the 
conversion  from  barotropie  energy  to  baroclinie  energy  The 
energy  conversion  favors  diurnal  tides  over  semidiurnal 
tides.  The  ridge  slope  controls  the  generation  of  linear 
internal  tides  from  a  ridge  but  energy  needed  to  produce 
large-amplitude  internal  tides  eomes  from  lee  waves.  If  the 
Froude  number  is  not  too  small  on  a  supercritical  ridge,  lee 
waves  enhance  the  westward  propagating  internal  tides  to 
form  tidal  bores  and  ISWs. 

[39]  The  second  step  of  ISW  generation  requires  a  sharp 
thcrmochnc  to  trap  energy  in  vertically  propagating  wave 
beams  to  form  horizontally  propagating  internal  tides.  When 
the  thcrmoclinc  is  diffuse  or  erased  by  the  presence  of  a 
deep  mixed  layer,  horizontally  propagating  internal  tides  arc 
weak,  and  no  ISWs  form.  A  shallow  upper  ocean  mixed 
layer  above  the  thcrmoclinc  does  not  affect  ISW  formation. 
The  mixed  layer  southwest  of  Taiwan  is  shallow  from  April 
to  August.  Not  surprisingly,  satellite  observations  reveal 
ubiquity  of  ISWs  during  this  period.  ISWs  arc  mostly  absent 
cast  of  the  Luzon  Strait,  where  a  sharp  upper  ocean 
thcrmoclinc  is  missing. 
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